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In certain Mott-insulating dimerized antiferromagnets, triplet excitations of the paramagnetic 
phase can decay into the two-particle continuum. When such a magnet undergoes a quantum phase 
transition into a magnetically ordered state, this coupling becomes part of the critical theory pro- 
vided that the lattice ordering wavevector is zero. One microscopic example is the staggered-dimer 
antiferromagnet on the square lattice, for which deviations from 0(3) universality have been re- 
ported in numerical studies. Using both symmetry arguments and microscopic calculations, we 
show that a non-trivial cubic term arises in the relevant order-parameter quantum field theory, and 
assess its consequences using a combination of analytical and numerical methods. We also present 
finite-temperature quantum Monte Carlo data for the staggered-dimer antiferromagnet which com- 
plement recently published results. The data can be consistently interpreted in terms of critical 
exponents identical to that of the standard 0(3) universality class, but with anomalously large cor- 
rections to scaling. We argue that the two-particle decay of critical triplons, although irrelevant in 
two spatial dimensions, is responsible for the leading corrections to scaling due to its small scaling 
dimension. 



I. INTRODUCTION 

Coupled-dimer magnets have become model systems 
for quantum phase transitions (QPT)i 2 . Such dimerized 
magnets are obtained from placing quantum spins on a 
regular lattice in d spatial dimensions, with two spins per 
unit cell and strong (weak) exchange interactions within 
(between) the unit cells. Depending on the ratio of the 
exchange interactions, the ground state can be either a 
paramagnct, dominated by singlet pairs in each unit cell, 
or a state with magnetic long-range order. Experimen- 
tally, the QPT between these two phases can often be 
driven by pressure. Additionally, these systems show an- 
other QPT upon applying an external field to the para- 
magnct , which allows to realize Bose-Einstein condensa- 
tion of magnonsi^ 

On the theoretical side, it is commonly assumed that 
the zero-field QPT is - by virtue of quantum-to-classical 
mapping - in the same universality class as that of 
the (d + l)-dimensional classical Heisenberg model, of- 
ten referred to as 0(3) universality class (note that the 
dynamical exponent z = 1 in that case)»^ Numerical 
simulations of various microscopic two-dimensional (2d) 
coupled-dimer Heisenberg models have indeed found crit- 
ical exponents consistent with three-dimensional (3d) 
0(3) universality, in agreement with this prediction^— 
Therefore it came as a surprise when results from accu- 
rate quantum Monte-Carlo (QMC) simulations of a par- 
ticular coupled-dimer Heisenberg magnet, the so-called 
staggered-dimer model with spins 1/2, displayed distinct 
deviations from standard 0(3) critical behavior, indicat- 
ing a different universality class ^ In contrast, the often 



studied columnar dimcr model was found to follow 0(3) 
universality^ suggesting the existence of two classes of 
coupled-dimer magnets* 1 ^ Subsequent QMC simulations 
of the staggered-dimer model^ 1 - - — focusing exclusively 
on the correlation-length exponent, obtained data con- 
sistent with those of Ref. [it However, it was argued that 
the data for the largest systems could be fitted to stan- 
dard 0(3) scaling laws. 

In this paper, we propose a resolution to the puzzle 
provided by the numerical data. We show that there 
are indeed two different classes of coupled-dimer mag- 
nets, henceforth called A and B. While class A follows 
standard 0(3) universality, the low-energy quantum field 
theory of class B is characterized by an additional cubic 
term which describes two-particle decay of critical fluctu- 
ations and has no classical analogue. While similar cubic 
terms have appeared before in the literature in different 
contexts? 1 ^— their effect on the critical behavior has not 
been discussed to our knowledge. We also note that the 
two-particle decay of triplet excitations at elevated ener- 
gies has received some attention both experimentall y 19 ' 20 
and theoretically) 2 1 ' 22 but it was commonly assumed that 
such processes are negligible at lowest energies. Here, we 
derive and analyze the critical two-particle decay term of 
class B in some detail. While its precise characterization 
in two space dimensions presents a challenge, our results 
are consistent with this term being weakly irrelevant in 
the renormalization-group (RG) sense. 

This leads us to suggest the following scenario for class- 
B coupled-dimer magnets like the staggered-dimer model: 
The asymptotic critical exponents are the ones of the 
0(3) universality class, but anomalously large correc- 
tions to scaling arise from the two-particle decay term. 
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This scenario is consistent with the numerical data re- 
ported in Refs. In addition, we also present finite- 
temperature QMC results for the temperature scaling of 
the quantum critical uniform susceptibility, which lend 
further support to this scenario. 



A. Overview of results 

The QPT between a non-symmetry-brcaking param- 
agnetic and a collinear antiferromagnetic phase in an in- 
sulating magnet with SU(2) symmetry is typically de- 
scribed by a quantum field theory of the </> 4 type, 
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in standard notation. Here, (p a (x,T) is a 3-component 
vector order-parameter field describing magnetic fluctu- 
ations near the ordering wavevector Q, with a = x,y, z. 
For simplicity, the action has been written for real (p a 
(appropriate for time-reversal invariant Q) and isotropic 
real space; the generalization to other cases is straight- 
forward. The critical behavior of the model ([T]) is known 
to be of standard (d + l)-dimensional 0(3) universality. 

We show that the following spatially anisotropic cubic 
termiiriS 



*7o / d d rdT<p • (d x if x d T tp) , 



(2) 



with x being a particular space direction, appears in the 
low-energy field theory for 2d coupled-dimer magnets be- 
longing to class B. This term causes two-particle decay of 
critical fluctuations; it bears some superficial similarity 
with Berry-phase and winding-number terms, to be dis- 
cussed below, however, its prefactor 70 is not quantized 
and the field <p is not restricted to unit length. 

Our detailed analysis suggests that the cubic term 
S3 in d = 2 space dimensions is irrelevant in the RG 
sense, albeit with a small scaling dimension. It consti- 
tutes the leading irrelevant operator at the critical fixed 
point. Consequently, the asymptotic critical behavior is 
of 0(3) type, but with anomalous corrections to scaling. 
We show that this scenario is consistent with the existing 
numerical data. 



B. Outline 

The body of the paper is organized as follows: In Sec.lHl 
we introduce the microscopic models under considera- 
tion, together with the bond-operator representation of 
their Hamiltonian. We discuss the conditions for the oc- 
currence of cubic terms in the microscopic bond-operator 
formulation. Together with the knowledge of the mag- 
netic ordering wavevector, this allows to sub-divide the 
models into classes A and B, where a cubic term docs 



(B) or does not (A) occur in the low-energy field the- 
ory. Sec. IIIII is devoted to a careful derivation of this 
low-energy field theory for the magnetic ordering tran- 
sition in the presence of cubic terms of a specific model 
belonging to class B. The critical behavior of this field 
theory will in turn be discussed in Sec. IIVI We employ 
both scaling arguments and direct classical Monte Carlo 
simulations to assess the relevance of the two-triplon de- 
cay term, with the conclusion that the most plausible 
scenario is its weak irrelevancy. This conclusion is sup- 
ported in Sec. |V] by QMC results obtained for various 
coupled-dimer models of classes A and B. A summary 
and outlook will close the paper. Various technical de- 
tails are relegated to the appendices. 



II. LATTICE MODELS OF DIMERIZED 

MAGNETS 

In this paper, we consider Heisenberg models with 
spins 1/2, Sj, placed on a regular lattice with spatially 
modulated couplings. The general Hamiltonian is thus 



H 



<jj') 



(3) 



where the sum is over all pairs of lattice sites jj'. Specifi- 
cally, all models have two sites per unit cell, with an anti- 
ferromagnetic intra-cell coupling J' defining the dimers. 
Spins in different unit cells are connected by couplings 
J according to the underlying lattice geometry. We re- 
strict our attention to lattices without geometric frus- 
tration, where the classical ground state is unique up to 
global spin rotations. Various 2d examples, namely the 
staggered and columnar dimer models as well as the her- 
ringbone and bilayer model, are shown in Fig. [1] 

Quite generically, these coupled-dimer models possess 
a paramagnetic ground state for J' 3> J, without sym- 
metry breaking of any kind and dominated by intra-cell 
singlets (Fig. [1]). In contrast, for J' J, a scmiclassi- 
cal Neel state with broken SU(2) symmetry is realized. 
The critical properties of the resulting QPT as function 
of J j J' are the subject of this paper. 



A. Bond-operator representation 

An efficient microscopic description of the excita- 
tions of coupled-dimer models is provided by the bond- 
operator representation^ Switching to a lattice of dimer 
sites i, the four states of a dimer i can be represented us- 
ing bosonic bond operators {s|,4 Q } i a = x jUi z )i which 
create the dimer states out of a fictitious vacuum. Ex- 



plicitly (and omitting the site index i), \s) 



- „t 



10), 



\a) = 4J0), where \s) = (\U) - |+t»/V2, \x) = 
(-Itt) +\U))/V2, \y) = i(|tt) +\U))/V2, \z) = 
(\U) + lit) )/V2. The Hilbert space dimension is con- 
served by imposing the constraint s\si + ^2 a t ia ti a = 1 
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as 




(d) 




FIG. 1: (Color online) Two-dimensional coupled-dimer mag- 
nets considered in this paper. In all panels, thick (thin) 
bonds refer to Heisenberg couplings J' (J), solid (open) cir- 
cles represent spins Sj (Sf ) of each dimer. In addition, 
the singlet configurations in the paramagnetic ground states 
realized for J' 3> J are shown, a) Staggered-dimer, b) 
columnar-dimer, c) Herringbone-dimer, and d) bilayer Heisen- 
berg model on the square lattice. The QPTs to the antiferro- 
magnetic phases are located at a) (J'/J) c = 2. 5196(2)i 
b) (J'/J) c = 1.9096(2)£^ c ) (J'/J) c = 2.4980(3)i2 d) 
(J'/J) c = 2.5220(1).— From the analysis in this paper, we 
conclude that the QPT of the models a) and c) belong to 
class B, while that of b) and d) belong to class A, for details 
see text. 



H = E + n 2 



H3 



T~Li + "Hrcst j 



(5) 



where Eq = 3J'N/8, with TV being the number of sites of 
the original (spin) lattice, and W.2,3,4 are terms obtained 
from ~H containing two, three, and four triplet operators, 
respectively. The approaches (a) and (b) only differ in the 
treatment of those higher-order terms in "H r cst which orig- 
inate from the different treatment of the Hilbert-space 
constraint: while (a) leads to an infinite series (starting 
at quartic order), in case (b) "H r cst only consists of the 
infinite on-site (i.e. hard-core) repulsion of triplets^ 

Performing a Fourier transformation for the triplet op- 
erators, ft = N'^ 1 / 2 yVexp(— ik ■ RAtt , here on the 
lattice of dimcrs with N' = N/2, can ^ c gencri- 

cally written as 



n 2 = V A r tt i,- + -B t 



tt tt +H.c. 

ka ka 
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on every site i. The original spin operators S 1 and S 2 of 
each dimer are given by 



SI 



2 



+ tl-s =f teapot-/) 



(4) 



where e a p-y is the antisymmetric tensor with e xyz = 1 and 
summation convention over repeated indices is implied. 

Using Eq. (j4]), the Heisenberg Hamiltonian ([3]) can 
now be re-written in terms of the bond operators 
{sj, ti a }2£-~— In the paramagnetic phase it is desirable 
to have a theory for triplet excitations only. Among oth- 
ers, two routes have proven useful: (a) In the spirit of 
spin-wave theory, the constraint is resolved by writing 
Si = s\ = (1 — tfatia) 1 / 2 . Expanding the root then gen- 
erates a series of higher-order triplet terms, (b) The 
formalism is re-interpreted as follows>2£ Starting from 
a background product state of singlets on all dimers, 
IV'o) = IL s l the operators t\ a can be viewed as cre- 
ating local triplet excitations in the singlet background. 
This re-intcrprctation e.g. changes the local triplet en- 
ergy from J'/4 to J 1 = J'/4 - (-3 J'/4). The constraint 

takes the form of a hard-core condition, Y^aAa^ia < 1- 
In both cases, subsequent approximations are usually de- 
signed to describe dilute triplet excitations on top of the 
paramagnetic ground stated 

The final bond-operator Hamiltonian can be written 



with the coefficients A^, Bj:, and 7g depending on the 
lattice geometry. Their explicit form for the models in 
Fig. QJ,,b will be given below. 

In the paramagnetic phase, an expansion around the 
singlet product state \?po} is justified. The leading-order 
term, r H,1^ describes non- interacting triplet excitations 
("triplons"), with energy 



t = J A 2 - B 2 , 



(9) 



obtained from a Bogolinbov transformation tt 

ka 



u k b l 

K k< 



k —ka 



with coefficients 
A n /2u> u ,-l/2 + A u /2uj n and u % v % = B^/2^. 
paramagnetic phase, ujj: > for all k. Setting H 
often referred to as harmonic approximation. 



1/2 + 
In the 
s %2 is 



H3, 4. rest contain interactions among the triplons. 
While H4 and "H r ost will contribute to the quartic self- 
interaction in the low-energy field theory, the cubic term 
H3 requires a more detailed discussion: as will be shown 
below, it may induce a cubic term of the form ([2]) in 
the low-energy theory, which then allows two-particle (in 
addition to three-particle) decay of critical fluctuations. 
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B. Bond operators for the staggered and columnar 
Heisenberg models 

For the staggcrcd-dimer model in Fig. QJ,, a straight- 
forward calculation gives 

H = J' + B % , 

B% = -— [cos(2fc x ) + cos{k x + k v ) + cos(k x - k y )] , 

Cj: = - J [sin(2fc x ) + sin(fc x + k y ) + sm(k x - k y )] , 

1% = B % ( 10 ) 

where k X:V refer to momenta on the original square lattice 
of spins. The triplon dispersion obtained at the harmonic 
level, i.e., from "H 2 > is shown in Fig. [2Jt. Its minimum 
energy at Q = (0,0) reaches zero at the critical value 
J' / J = 3; the quantum Monte Carlo result for the loca- 
tion of the QPT is {J'/J) c = 2.5196(2)£± 2 . 

Similarly, for the columnar dimcr model, Fig. [Tp, 

H = J' + B % , 

Bj: = — [2cosfc. y - cos(2k x )} , 

= -J sm(2k x ), 

7j: = — g [2 cos k y + cos(2fca;)] . (11) 

Now the dispersion minimum is at Q = (0, ±7r), Fig. [2p. 
At the harmonic level, the critical point is again located 
at J' I 'J = 3, while the currently most precise quantum 
Monte Carlo result is (J'/J) c = 1.90948(4) ,22, consistent 
with the previous value of (J'/J) c = 1.9096(2)4 

Fig. [2] also shows the lower bound of the two-particle 
continuum at the critical coupling. In the columnar 
dimer model, Fig. [5p, two-particle decay is only possible 
at elevated energies (where the single-particle dispersion 
is inside the two-particle continuum). In contrast, in the 
staggered-dimer model the single-particle dispersion co- 
incides with the bottom of the two-particle continuum at 
criticality. While this coincidence at all wavevectors is 
an artifact of the harmonic approximation, it is the cor- 
rect result near Q = for a critical point with dynami- 
cal exponent z = 1. In the presence of a non- vanishing 
cubic term H3 in the Hamiltonian, this implies a cou- 
pling between one-particle and two-particle sectors down 
to lowest energies, as discussed in the following. 

C. Symmetries and two-particle decay 

Before diving into the derivation of the low-energy held 
theory, it is worth discussing for which coupled-dimer 
models a cubic term, describing two-particle decay of 
triplons, will be part of the low-energy field theory. 

First, a cubic piece "H 3 with £g 7^ does not exist in 
all microscopic models. In fact, in high-symmetry cases 
like the much-studied bilayer Heisenberg model, Fig. [T}1, 



£g = 0i2I An analysis shows that, due to the antisym- 
metric character of H3, its coefficient £g is non- vanishing 
provided that two dimers i, i' are coupled in an asym- 
metric fashion, such that the couplings J kk (fc, k' = 1, 2) 
between the spins S k on i and S k on %' obey J 11 — J 22 7^ 
and/or J 12 — J 21 7^ 0. This can be translated into the 
following symmetry condition: £g vanishes provided that 
the model remains invariant if in every dimer the spins 
1 and 2 are inter-changed (together with all their cou- 
plings). 

Second, H3 can enter the low-energy theory only if the 
ordering wavevector on the dimer lattice is Q — - this 
immediately follows from momentum conservation. As 
we shall show in the next section, these two conditions 
are indeed sufficient for a non-vanishing cubic term J2]) 
to appear in the low-energy theory, and thus these con- 
ditions define our class B of coupled-dimer models. 

Members of class B are the staggered-dimer model in 
Fig. [Th as well as the coupled-dimer models on the hon- 
eycomb lattice^ and the herringbone square lattice^ 
Fig. [T}; - these are exactly the models for which devia- 
tions from 0(3) critical behavior have been discussed 
On the other hand, both the columnar and bilayer mod- 
els, Figs, [lb and d, do not fulfill these conditions and 
hence belong to class A. Indeed, standard 0(3) critical 
behavior has been established for these modelsi£&2i 



III. LOW-ENERGY QUANTUM FIELD 
THEORY 

In this section we derive the effective low-energy field 
theory designed to capture the physics near the antiferro- 
magnetic quantum critical point of the models introduced 
in Sec. [XXI We shall present in detail the derivation of a 
ip 4 theory from the bond-operator representation of the 
microscopic spin-1/2 model; in Appendix[B]we shall also 
sketch the derivation of a non-linear sigma model in the 
scmiclassical limit of ([3]). In both cases, a cubic term of 
the form will appear for the class-B models; however, 
we believe the non-linear sigma model is not useful for 
further analysis, see Sec. MI Bl 

A. From bond operators to the tp 4 model 

A derivation of a <p theory from the bond-operator for- 
malism has been presented in the context of the columnar 
dimer model in Ref. [29|. Here we shall follow this proce- 
dure, taking also into account the cubic piece H3 which 
will enter the ip 4 theory only if the ordering wavevector 
Q = 0. 

In a Lagrangian formulation, the bond operators can 
be represented by a complex bosonic vector field i(r, r). 
Importantly, this contains information both about the 
staggered and uniform magnetization fluctuations on 
each dimer, i.e., it contains the degrees of freedom of 
both (S 1 - S 2 ) and (S 1 + S 2 ). The latter live at high 
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FIG. 2: (Color online) Triplet excitation spectrum (harmonic level) for the staggered (a) and columnar (b) dimerized AF 
Heisenberg models along high symmetric directions in the Brillouin zone. Solid (black) and dashed lines correspond respectively 
to J' = 3J and J' = 3.5J. The dotted line is the bottom of the two-particle continuum for J' = 3J - this coincides with the 
single-triplet dispersion in (a). Finally, panels (c,d) show the Brillouin zones for the staggered (c) and columnar (d) dimerized 
Heisenberg models. 



energies and will eventually be integrated out to obtain 
a theory for the staggered fluctuations only. To this end, 
we decompose the complex field t = Z{tp + iaif), where ip 
and if are real three-component vectors, a is the lattice 
spacing, and 



Z = 2~ 



• i J- 



(12) 



is a renormalization factor. 

A continuum-limit formulation is obtained by an ex- 
pansion of momenta in the vicinity of k = Q. After 
some straightforward algebra the action corresponding 
to Hamiltonian W 2 © takes the form 

1 



S 2 = -jd 2 rdr[cl{d x 0) 2 +c 2 y {d y p) 2 +m o p 2 ] 

+ — J d 2 rdr {m^if 2 + 2i s Jm lx Tfd T (p) (13) 



for both the staggered and columnar dimer models. Here, 
mo tunes the phase transition, while m^ remains finite 
near the QPT. Explicitly, we have mo = J' (J' — 3 J) a~ 2 
and m-H = J 12 . This implies that the transition at the 
mean-field level takes place at J' = 3J (as in the har- 
monic bond-operator approximation above). From the 
form of mo we can deduce its bare scaling dimension to 
be 2 (relevant), since the continuum limit of the field 
theory obtains for a — > 0, and the mass term grows upon 



approaching that limit. Furthermore, 



3JJ' and 



Cy = J J'. The dynamics is encoded in the mixed- field 
linear time derivative. 

The cubic Hamiltonian piece ([7]) of the staggered- 
dimer model, Fig. QJi, can be cast into the form 



S3 = 7o / d 2 rdrd x ip ■ {(p x if) 



(14) 



where d x originates from £g oc k x for small k, and 70 = 
JJ'f 25ct3. The latter implies that the scaling dimension 



of 70 is (—1/2), and consequently it is irrelevant at tree 
level. Note that there is another term of the form 



S3 



7 / d rdr d x if ■ {ip x 7?) 



(15) 



This term turns out to be more irrelevant than £3 (it 
scales like a 5 ) and hence will be discarded. 

Importantly, a term of form £3 (|14p does not appear for 
the columnar dimer model, Fig. QJ>: there, Q = (0,±7r), 
and a term with three fields carrying momentum Q is 
forbidden by momentum conservation. In contrast, for 
other coupled-dimer models with ordering wavevector 
Q = and non- vanishing in ^3, a cubic term ap- 
pears. We have checked this explicitly for the honeycomb 
and herringbone lattices. The spatial structure of 1S3 is 
strongly anisotropic, with the direction of the scalar spa- 
tial derivative being determined by the orientation of the 
dimers (via the small- momentum expansion of For 
instance, for the herringbone lattice, Fig. [TJ;, this direc- 
tion is diagonal w.r.t. the underlying square lattice. 

The quartic interaction term is again identical for the 
staggered and columnar dimer models and reads 



S. 



^ / d 2 rdr 
24 



{ip^Y + ai (<p x jf) + a 2 (p 2 TT 2 ,(16) 



where uo = 3 • 2 d+3 J' 3 a~ 1 and ai,a 2 oc a 2 . Since the 
terms cx a 1,02 are more irrelevant than the first one we 
discard them. 

The next step is to integrate out the field if. Although 
jtItt > 0, care is required: the coupling to (p via S3 renders 
the field if gapless at the critical point. However, we 
have explicitly checked that this complication does not 
introduce singularities, see Appendix [X] To proceed, we 
introduce a new field if ' = if + —j=d T p, such that the 
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quadratic part of the action reads 

S 2 = \ I SrdT [c 2 x (d x 0) 2 + c 2 y (d y 0) 2 + (0 T p) 2 + m o 2 ] 



(17) 



where now a second-order time derivative for ip appears. 
Next, 7?' is integrated out from S2 + S3 + S4. Keeping 
only the lowest-order terms and defining £24 = S2 + S4 
we find 

= I I d 2 rdr [cl(d x 0f + c 2 (d y (p) 2 + (d T 0) 2 + m o¥ ? 21 



^jd 2 rdr(0 2 y 



(18) 



S3 = ijo / d 2 rdr ip ■ (d x x d T (p) . (19) 



and 



The additional term S3 p^|) - identical to announced 
in the introduction - represents the crucial difference 
to a standard ip 4 (or Ginzburg-Landau) theory as de- 
scribed by 1S24. Among the various additional higher- 
order terms, there are 



S3+2n oc i d 2 rdr p ■ (d x p x d T ip) ip 



■2ti 



(20) 



with n > 1. As usual, they may be discarded since they 
are more irrelevant in the RG sense compared to £3 (flU|) . 

Summarizing, within the framework of the bond- 
operator approach one can derive, for the staggered dimcr 
model, an effective ip 4 theory supplemented by an infinite 
number of additional terms. The most relevant of these 
is S3 (fT9"| which will be further discussed in Sec. IIII Bl 
below. 



B. Discussion of the cubic term 

The derivation of the ip 4 theory has lead to a cubic 
term £3, Eq. (|T9)) . which is present for the staggered 
dimer model, but absent for the columnar dimcr model. 



1. Symmetries 

The cubic term (fH?)) respects SU(2) symmetry and time 
reversal, as is odd under time reversal. It respects 
momentum conservation provided that (p parameterizes 
fluctuations at wavevector zero. Finally, the existence 
of the term requires that (p is odd under the mirror op- 
eration x — > — x while it is even under y — > — y, which 
is the case for the staggered magnetization on the hori- 
zontally aligned dimers in the models in Fig. [TJi,b. The 
cubic term is fundamentally quantum in the sense that 
no quantum-to-classical mapping exists for this term due 
to its prefactor i. (This also implies that the field theory 
with cubic term is not amenable to an efficient Monte 
Carlo sampling, since it suffers from a sign problem.) 



2. Quantization and relation to topological charge 

A cubic term can also be derived in the language of the 
non-linear sigma model, see Appendix [B] with the result 

S3 = iSo / d d rdr n ■ (d x n x d T n), (21) 



where n is now a unit-length 0(3) field, and 6 is pro- 
portional to the modulation (J — J') of the couplings in 
Fig. QJl Similar cubic terms were derived before for the 
standar d 14 i 15 and the dimerizecUi Heisenberg models on 
the honeycomb lattice, but in all cases neglected in the 
subsequent analysis. 

Importantly, Eq. (f2"Tj) may suggest an interpretation 
in terms of a topological charge (or skyrmion number) in 
x-t space. We can introduce a functional Q of a vector 
field a in two dimensions according to 



Q[a(x,y)] = J dxdya- (d x a x d y a). 



(22) 



For a unit-length field a, Q[a\ is known as topological 
term. It is quantized to integer values for periodic 
boundary conditions and smooth configurations of a: a 
is a map from the 2d plane to the unit sphere, and Q 
measures how often space is wrapped around the sphere.— 
Notably, the Berry phase term Sb of a Id antiferromag- 
netic spin chain, represented in spin coherent states, can 
be represented in a very similar fashionj^^ 



1— j dxdrn- (d x nx d T ft) = 2iriSQ[n(x,T)]. (23) 



Due to the quantization, Sb drops out from the partition 
function for integer spins S due to e SB = 1, while Sb 
contributes non-trivial sign changes from skyrmions for 
half-integer Now, the dimer-model cubic term (|2"Tj) 
in d = 2 can be written as: 



S 3 = 4iri5 / dy Q[n y (x,T)} 



(24) 



where ft y (x, r) = n(x, y, r). Note that, in contrast to Sb 
in (|23p . the prefactor in £3 ([24j) is not quantized. Based 
on the expression (|24|) , Ref. [l7| concluded that S3 is neg- 
ligible, arguing that smooth configurations imply that 



Q[f 



const, and skyrmion lines described by non-zero 



Q[n y ] should be energetically suppressed. This argument 
is certainly not rigorous, as instanton events are not ac- 
counted for. However, the non-universal prefactor of S3 
may suggest that non-trivial contributions to Q[wj,] tend 
to average out. 

A central issue in the discussion of S3 is therefore 
whether the quantization in terms of a topological charge 
in x-t space really plays a role. We believe that this 
is not the case, for the following reasons: (i) In the ip 4 
(i.e. soft-spin) version of the field theory, the field oc- 
curring in S3 (|19p is not normalized to unity, such that 
Q^y] is not quantized. (Even if amplitude fluctuations 
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are frozen out at large length scales, Q[0y) is sensitive to 
fluctuations on all scales.) (ii) Before taking the spatial 
continuum limit, the expression in 1S3 involves a discrete 
sum over x, with the derivative d x replaced by a linear 
function of 0. While lattice definitions of Q preserving 
its topological character for unit-length fields have been 
put forward^ those involve the fields in a strongly non- 
linear fashion. In contrast, for a discretization with a 
linear approximation to the derivative it is easy to show 
that the topological character is not preserved. This will 
be explicitly shown in Sec. IIVBI below. 

Therefore, we believe that the cubic term in the order- 
parameter theory of class-B couplcd-dimcr magnets is un- 
related to quantized topological charges, i.e., the relation 
suggested by Eqs. (f2~Tj), (|24|) is an artifact of the unit- 
length continuum limit underlying the non-linear sigma 
model. Consequently, standard tools like perturbative 
RG can be used to analyze the cubic term in the ip 4 for- 
mulation. 



IV. CRITIC ALITY IN THE PRESENCE OF 
TWO-TRIPLON DECAY 

According to our analysis so far, the cubic term £3 
is the most relevant additional term present in the low- 
energy field theory for the quantum phase transition (as 
compared to the standard 0(3) case). Therefore, the cen- 
tral question is whether this term is relevant or irrelevant 
in the RG sense at the 0(3) (or Wilson-Fisher) critical 
fixed point in (2 + 1) dimensions. This can be answered 
by determining the scaling dimension of the coupling con- 
stant multiplying the local cubic operator which appears 
inside £3. 

In this section, we shall follow two routes: First, we an- 
alyze the cubic operator in the <p theory perturbatively 
in e = 4— D. Second, we determine the operator's scaling 
dimension directly in D = 3 dimensions, by means of a 
classical Monte Carlo calculation of the operator's cor- 
relation functions at the Wilson-Fisher fixed point. The 
results of both methods are consistent with the cubic 
term being weakly irrelevant in D = 3. 



Perturbative determination of the scaling 
dimension 



The 0(3) critical fixed point is perturbatively acces- 
sible in the framework of the </? 4 theory ([1]) in a double 
expansion in the quartic interaction uq and e = 4 — D. 
For D < 4 the Gaussian fixed point is unstable towards 
the Wilson-Fisher fixed point, with a renormalized inter- 
action u ~ 0(e). 

We start by determining the scaling dimension of the 
cubic operator's ([2]) coupling constant at the Gaussian 
fixed point. After re-scaling the lengths such that the 
gradient terms are isotropic, the action in D = d + 1 



dimensions reads^ 



5 



m ip a 2 + {V(p a y 



+ i7o / d D r 0- {d x x d y 0) 



U ° / AD (2^ 

(25) 



At tree level, we obtain the well-known scaling dimen- 
sions 

[0] G = (D-2)/2, 
[u ] G = D-4[0] G =4-D, 
[ 7o ] G = D - 2 - 3[<^] G = (2 - D)/2, (26) 

with the subscript G referring to the Gaussian fixed 
point. Substituting D = 3, the cubic term is found to 
be irrelevant with a scaling dimension of [70] = — \ - the 
same conclusion appeared already in Sec. IIII Al 

At the Wilson-Fisher fixed point, both fields and ver- 
tices receive perturbative corrections leading to anoma- 
lous dimensions. A simple (but incomplete) estimate of 
the scaling dimension of 70 at the Wilson-Fisher fixed 
point consists of taking into account the field renormal- 
ization only. This amounts to using [70] = D — 2 — 3[0] 
with [0] = (D - 2 + rj) /2 leading to 



[7o] 



D 



377 
2 



-0.55625 



(27) 



where 77 = 0.0375(5) in D = 3 (Ref. 133J) was used. Al- 
though indicative, we cannot expect this estimate to be 
reliable, as it ignores vertex corrections: it is known that 
composite operators may have large anomalous dimen- 
sions (see, e.g., Ref. [13). 

A more complete treatment requires a perturbative RG 
analysis of the full theory S24+S3. This expansion is done 
about the Gaussian theory, with two dimensionless non- 
linear couplings u = uqK d ~ a and 7 = 7 oA — ~ , where 
A is an ultra-violet cutoff. To one-loop order, the calcu- 
lation is conveniently performed in the momentum-shell 
scheme. It turns out that, due to the antisymmetry of the 
7 vertex, no diagrams mixing u and 7 exist to one-loop 
order. Furthermore, 7 does not introduce field renormal- 
izations. Hence, the flow equation for u is not modified 
by 7, and the flow of 7 does not involve u. To one-loop 
order we simply have: 



du 



(A-D)u-K d ^±^u\ 



2 - D 



(28) 
(29) 



where dl = rfA/A, N = 3 is the number of field com- 
ponents, and K d = (2 d - 1 7r d / 2 T(d/2))- 1 . Thus, the tree- 
level result [70] = (2 — D)/2 does not receive one-loop cor- 
rections. If rcnormalizations of the 7 vertex due to u re- 
mained absent at higher loop orders, only field rcnormal- 
izations would influence the flow of 7, and the estimate 
(|2T|) would be correct. However, we see no fundamental 
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reason for a general cancellation of such vertex renormal- 
izations. Instead of going to higher loop orders, we will 
improve on the estimate (|27j) using a non-pcrturbative 
numerical approach. 



B. Monte Carlo analysis in D = 3 

We shall now numerically determine the scaling dimen- 
sion of the cubic term 1S3 directly in (2 + 1) dimensions 
at the Wilson-Fisher fixed point. Note that this task is 
simpler than solving the full quantum model including 
53: in particular, it boils down to the simulation of a 
classical problem in D = (d + z) dimensions, with z = 1, 
as the 0(3) critical field theory described by S24 follows 
a quantum-to-classical mapping. 

We define the composite operator 



0(f) = tp{r) ■ (d x tp(r) x d y tp{r)) 



(30) 



Its scaling dimension [0] = A<p can be obtained from the 
long-distance decay of its correlation function: 



C*(f) = (0(f)0(O)) <x 



1 




12A f 



(31) 



FIG. 3: Histograms of the intra-layer "Skyrmion number" 
Q(z) (see text), obtained from a Monte-Carlo simulation of 
the classical Heisenberg model (|33[) at its critical tempera- 
ture. The curves have been obtained from 10 8 measurements 
on systems of size L 3 with L = 4,6, 8; each measurement 
gave Q(z) for a single layer with fixed z. Note that Q is not 
quantized, as it does not involve unit-length fields in the con- 
tinuum limit. The distributions are found to be Gaussian, 
with a width scaling linearly with L. 



From this, the scaling dimension of the coupling con- 
stant (more correctly: the associated vertex function) is 
obtained through 



[70] = D - A 



o 



(32) 



In the following, we determine the scaling dimension 
Aq of the composite operator by a lattice Monte 
Carlo simulation of a classical Heisenberg ferromagnet 
in D = 3 dimensions, where we shall measure the corre- 
lator Eq. (f3~Tj) at criticality. This approach exploits that 
the model is in the same universality class as the 0(3) 
Landau-Ginzburg theory and hence realizes the Wilson- 
Fisher fixed point in D — 3, but gives us access to correla- 
tion functions in a non-perturbative manner. Specifically, 
we simulate the classical Heisenberg model 



H — — J S Si ■ Sj 



(33) 



<ij> 



with ferromagnetic interactions between nearest neigh- 
bors on a simple cubic lattice. The Si are classical (com- 
muting) three-component vectors of unit length (Sf = 
1). We employ the Wolff cluster algorithm,— which al- 
lows an efficient Monte-Carlo simulation and provides 
high-accuracy critical exponents for the 0(3) universal- 
ity class^ The critical point of this model is known to 
be located at K c = J/(k B T c ) = 0.693035(37)^ 

In the lattice simulation, the operator needs to be 
discretized. Guided by the derivation of the field theory 
from the discrete lattice model, Sec. UH we know that 
the derivatives in Eq. ([3U1) should be discretized using a 



linear function of the spins (in contrast to Ref. |32| ) 
standard two-point forward formula leads to 



O 



i, lattice 



= Si ■ s 



x S 



i+e. 



The 



(34) 



where e x (e y ) denotes a unit step in x (y) direction. We 
have checked that other (linear) discretization schemes 
give qualitatively similar results. 

Before we turn to the results for the correlator C(r), 
we make a brief detour to discuss the quantity it- 
self. As mentioned above, its layer integral Q(z) = 
^2, X yO(x,y, z) / {At:) may suggest an interpretation in 
terms of a topological charge. However, the numerical 
Monte Carlo simulations show that Q(z) is not quan- 
tized at criticality, see Fig. [31 Instead, Q(z) displays a 
single peak at Q — 0, with a width scaling as L in a sys- 
tem of size L 3 . Since the number of spins in each layer 
is L 2 , this width simply reflects the standard thermody- 
namic scaling of fluctuations of a non-critical extensive 
observable. Therefore, Fig. [3] supports the conclusion of 
Sec. lIII Bl that 0(f) is a conventional non-critical density. 

The correlator of is measured along the two inequiv- 
alent directions, i.e., within the xy-plane and along the 



C xy (r) = (0(^0,0)0(0)) = (0(0, r, 0)0(0)), 
C z (r) = (0(O,O,r)0(O)) 



(35) 



where r now denotes discrete lattice coordinates. We 
find that both correlation functions drop quickly with the 
separation r, and a large number of Monte Carlo sweeps 
are required to reduce the statistical uncertainty in the 
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FIG. 4: The correlators C xy {L/2),C z (L/2) as functions of 
L/2 for different system sizes L from Monte-Carlo simula- 
tions of the ferromagnetic Heisenberg model (|33p at its crit- 
ical temperature. The dashed line corresponds to a decay 
proportional to (L/2)~ 6 . 



correlator. The most efficient way to estimate the decay 
exponents of C{r) in a finite-size system is to measure 
the correlation functions at half the linear system size, 
C xy (L/2) and C z (L/2), for different lattice sizes, as to 
minimize finite size effects. We employed up to 10 12 Wolff 
cluster updates for system sizes L = 6,8, 10, 12, 14 to 
obtain the data shown in Fig. 

Despite the relatively small system sizes, C xy (L/2) and 
C z (L/2) show an algebraic decay with L/2, consistent 
with critical behavior. For both correlators, the decay ap- 
pears to be faster than 1/r 6 . However, a reliable determi- 
nation of the decay exponent is difficult, because the data 
display a slight curvature at the largest system sizes. A 
direct fit of all data points yields Aq ~ 4, while the large- 
system data is more consistent with A<p rj 3. 2... 3. 5. 
With Eq. (|32|) this suggests that the scaling dimension 
[7o]w_f is in the range —0.2 . . . — 0.5, consistent with weak 
irrelevancy of the cubic operator at the Wilson-Fisher 
fixed point. However, from the present data we cannot 
rule out that 70 is instead weakly relevant. 



C. Scenario: Large corrections to 0(3) scaling 

Let us summarize the state of affairs concerning the 
critical behavior of class-B dimer models: (i) We have 
found that the combination of low symmetry and vanish- 
ing ordering wavevector leads to the presence of a cubic 
term in the low-energy field theory - this cubic term rep- 
resents the most relevant difference to a standard 0(3) 
field theory, (ii) The cubic term is strongly irrelevant in 
(3 — e) space dimensions, while the Monte-Carlo results 
of Sec. IIV Bl suggest that it has a small scaling dimension 
in d = 2, most likely being weakly irrelevant, (iii) The 
published QMC results for the staggered dimer model 
indicate cither critical exponents slightly different from 
those of the 0(3) universality class^ or, if the fitting is re- 



stricted to large systems only, exponents consistent with 
their 0(3) values! 11 ' 12 Note that a conventional value for 
the critical exponent v was obtained also using an uncon- 
ventional finite-size scaling analysis of the spin stiffness^ 
whose validity remains to be verified. 

Points (i) and (ii) strongly suggest that the cubic term 
is responsible for the unusual behavior seen in the QMC 
calculations. Point (iii) then implies that the cubic term 
is irrelevant (instead of relevant) in the RG sense, as 
otherwise the deviations from 0(3) universality would 
grow (instead of shrink) with system size. 

This leads us to propose the following scenario for the 
quantum phase transition in class-B coupled-dimer mod- 
els: The cubic term is weakly irrelevant in d = 2 and 
therefore constitutes the leading irrelevant operator at 
the Wilson-Fisher fixed point. Hence, the asymptotic 
critical behavior is that of the standard 0(3) univer- 
sality class, but the corrections to scaling are different 
from standard 0(3) universality. In the next section, we 
present numerical results that support this scenario. A 
detailed analytical calculation of the corrections to scal- 
ing arising from the cubic term is left for future work. 

V. QUANTUM MONTE CARLO RESULTS 

In this section, wc present further results from QMC 
simulations of critical dimerized antiferromagnets to as- 
sess the scenario discussed above. First, we re-analyze 
finite-size data for critical exponents obtained in Ref. 
in terms of anomalously large corrections to 0(3) scaling. 
Second, we show results for the finite-temperature uni- 
form susceptibility at criticality - this is a non-critical 
quantity which has been studied both analytically and 
numerically for critical antiferromagnets in the past. 
Again, dimer models of classes A and B are found to 
display distinctly different behavior. 

A. Correction exponent 

Following the idea of large corrections to scaling out- 
lined above, we investigate here whether the inclusion of 
appropriate scaling corrections yields critical exponents 
compatible with the 0(3) universality class for the stag- 
gered dimer model. In particular, we focus on the order 
parameter (i.e. the staggered magnetization) m s , esti- 
mated using only the z-component of the spin operator 
via m z s = I J2i(~^) x+v Si\- The analysis of this quantity 
lead to the most pronounced deviation from 0(3) critical 
behavior in Ref. [9J, and it complements more recent sim- 
ulation a 11 ' 12 that focus solely on the exponent v of the 
correlation length. At the quantum critical point, m z s is 
expected to scale according to 

(m;)~i-^(l + c m i- w ), (36) 

thus providing access to the ratio fi/v of critical expo- 
nents. For the case of the ladder model, perfect agree- 
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TABLE I: Fit results for the critical exponent quotient fi/u 
for the staggered dimer model. The table summarizes several 
results of fits including a correction to scaling exponent for 
three values of a c within twice its error bar and L > 10. The 
entry (ref) refers to the relevant reference values for the 3D 
0(3) universality class given in the last line. Reported error 
bars are twice the fit error. 



etc 






x7a\o.f 


2.5194 


0.525(1) 


ref 


4.0 


2.5196 


0.529(1) 


ref 


1.33 


2.5198 


0.533(1) 


ref 


1.1 


2.5194 


ref 


0.63(3) 


1.7 


2.5196 


ref 


0.55(2) 


0.6 


2.5198 


ref 


0.48(3) 


2.6 


Ref. [37 38 


0.5188(3) 


0.782(13) 





merit with 0(3) exponents was found — even when ne- 
glecting the presence of the corrections to scaling^^ 

For the standard 0(3) universality class, the correction 
exponent uj is given by21 ojq(3) = 0.782(13) (in the (p 4 
language arising from the flow of the quartic interaction 
u on the critical manifold). Hence, if the deviations in 
P/v observed in Ref. were due to standard irrelevant 
operators, inclusion of u = 0Jo(3) should result in the 
0(3) value for (3/v. Our results from performing fits to 
the data of Ref. [H, presented in Table Q] indicate that 
this is not the case, i.e., we cannot cast the fitting results 
with 0(3) values of the critical exponents. 

We continue to investigate a second scenario, in which 
we fix the known 0(3) exponent ft/v but leave w as a 
free fit parameter. The lower parts of Table U contain 
the corresponding fitting results, which indicate that we 
can indeed arrive at a 0(3) critical value for fi/v : but 
at the expense of a u < 0Jo(3)- It should be under- 
stood that the performed analysis actually provides an 
effective correction exponents for the length scales stud- 
ied. Nevertheless, from further numerical studies of the 
herringbone and honeycomb couplcd-dimcr models^ we 
can empirically relate the presence of the cubic operator 
to unusually large (and slowly vanishing) corrections to 
the leading 0(3) scaling. In fact, from the symmetry ar- 
guments presented in Sec. Ill C[ both models belong to 
class B (as the staggered-dimer model) with a cubic two- 
triplon decay term at low energies. This provides numer- 
ical support for the scenario outlined above, in which the 
cubic operator leads to enhanced corrections to an 0(3) 
scaling behavior. 



B. Finite-temperature susceptibility 

In this section, we focus on the thermodynamic behav- 
ior, in particular the temperature scaling of the uniform 
susceptibility at the quantum critical point. For this pur- 
pose, we performed QMC simulations of various dimcr- 
ized two-dimensional antiferromagnets at their respective 
quantum critical coupling ratios, see Fig. [T] 



For the simulations, we employed the stochastic se- 
ries expansion method with a directed operator-loop 
update^r— We considered systems with TV = 2L 2 spin, 
for linear system sizes L up to 512. The quantum non- 
linear sigma model prediction for the uniform suscepti- 
bility is a linear dependence x ~ AT on the tempera- 
ture (T) within the quantum critical region, where the 
prefactor A depends on the spin-wave velocity and a uni- 
versal constant i 42 ' 43 This implies an essentially constant 
ratio x/T ~ A inside the quantum critical region. Such 
a linear-T scaling of \ nas been observed for both the 
bilayer— and a coupled plaquette lattice 4 - model signifi- 
cantly into the quantum critical region. 

The QMC results of the uniform susceptibility for 
the staggered and the columnar dimer arrangements are 
shown in Fig. [5] for different system sizes. As seen 
from comparing the QMC results for different system 
sizes, we obtain finite-size converged estimates for the 
thermodynamic-limit behavior down to T/J = 0.04. 
Comparing the data for the two cases, we find that 
while for the columnar arrangement x/T shows only mild 
changes with T below about 0.2 J, for the staggered case, 
larger deviations from a constant value of x/T are ob- 
served over the whole accessible temperature range. We 
take these enhanced deviations from the lincar-T scaling 
of x as a signature of the scenario outlined in the pre- 
vious section, even though we are not in a position to 
derive from our theoretical analysis the actual form of 
the leading deviation from the linear-T scaling of X- 

In Fig. [6l we show the low-T temperature dependence 
of x/T for all four models in Fig. [T] at their respective 
quantum critical points, extrapolated to the thermody- 
namic limit. Strikingly, we find sizeable and similar cor- 
rections to the linear-T scaling of x for both class-B mod- 
els (staggered dimer and herringbone) , whereas such cor- 
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FIG. 5: QMC results for the temperature dependence of the 
uniform susceptibility \ for the columnar and the staggered 
model at the quantum critical point for different system sizes. 
Also shown are fits to an ansatz x/T = A — BT + CT 2 for the 
low-T behavior of x/T extrapolated to the thermodynamic 
limit (dashed). Note the logarithmic T scale. 
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FIG. 6: Quantum Monte Carlo results for the temperature 
dependence of the uniform susceptibility \ for different dimer 
antiferromagnets. The shown data represents the thermody- 
namic limit. Linear lines represent fits of the low-T behavior 
of x/T to a linear ansatz. 



rcctions are much less pronounced for the class-A models 
(columnar dimer and bilayer). 

For all models, the leading low-T behavior is consistent 
with a linear decrease of x/T, i.e., 

X/T = A — BT, B > 0, (37) 

shown by the linear fits in Fig. [5] For the staggered 
dimer and the herringbone model, enhanced corrections 
B are required, as seen from the linear fit lines in Fig. |6l 
which have about twice the slope as those for the other 
two models. While polynomial corrections to the linear- 
T scaling of x are thus compatible with our data, it is 
interesting to assess, if our data is consistent also with 
other functional forms of the leading correction terms. 
For example, recently Sandvik observed a dominating 
logarithmic term in the low-T corrections to the linear-T 
scaling of x in a two-dimensional Heisenberg model with 
four-spin interactions (the J — Q model) j 23 i 45 For our 
data, the type of corrections can be judged from both 
Figs. [5] and [6] Fig. O with a logarithmic T scale, shows 
a quadratic fit x/T = A — BT + CT 2 : which is seen to 
fit the data well up to about T w 0.15, while no robust 
behavior linear in log(T) is observed. 

We thus conclude from our QMC analysis, that (i) 
at the quantum critical point deviations from a linear- 
T scaling of x ar c exhibited by all considered models, (ii) 
the leading corrections to the linear-T scaling can be cap- 
tured by a low-order expansion in T fEq. l37[) . and (iii) in 
those models, for which non-trivial cubic terms emerge, 
considerably enhanced corrections to the linear-T scaling 
of x are present. In fact, Figs. [5] and [5] give a rather clear 
indication of the two classes A and B of dimer models, 
with class A (B) displaying small (large) corrections to 
the leading x/T = const behavior. It proved difficult to 
extract the actual functional form of the scaling correc- 
tions from the QMC simulations, but we tend to exclude 



low-temperature logarithmic corrections as found for the 
J-Q model in Refs. HUM 

VI. CONCLUSIONS 

Analyzing quantum phase transitions in models of 
couplcd-dimer magnets, we have identified two distinct 
universality classes A and B. While class A displays con- 
ventional 0(3) critical behavior, class B is characterized 
by the possibility of two-particle decay of critical fluctu- 
ations, described by a cubic term in the order-parameter 
field theory. We have shown that various 2d coupled- 
dimer models including the recently studied staggered- 
dimer model belong to class B. Combining field-theoretic 
arguments and results from large-scale numerical simula- 
tions, we have put forward the following scenario for the 
quantum phase transition in class-B models: The leading 
critical behavior is that of the standard 0(3) universal- 
ity class, but anomalously large corrections to scaling, 
different from 0(3) behavior, arise from the cubic term. 
This scenario appears consistent with all available infor- 
mation, in particular it solves the puzzle concerning the 
interpretation of recent QMC results for the staggered 
dimer modeller— A precise analytical characterization of 
the scaling corrections arising from the cubic term is left 
for future work. 

It is conceivable that similar two-particle decay terms 
also appear for quantum phase transitions with under- 
lying symmetries different from SU(2). Then, the cor- 
responding class-B transitions might even display novel 
leading critical behavior. 
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Appendix A: Derivation of if 4 theory: self-energy 
corrections 

In Sec. IIII A[ we derived a i^ 4 -type low-energy theory 
from the bond-operator representation of the staggered 
dimer Heisenberg model. In the course of the deriva- 
tion we integrated out the field tt based on the fact that 
it is gapped. However, at the critical point this state- 
ment turns out to be true only at tree level, due to the 
self-energy corrections to 7? arising from two-particle de- 
cay described by Eq. (H3J). Therefore we have to verify 
whether integrating out 7? is still permissible without in- 
troducing singular terms. 
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Consider the sclf-cncrgy of the field 7? due to the cubic 
term (| 14|) which enables decay into two <p fields. In lowest 
order perturbation theory and at criticality, it reads 



2 , 7T 

S ff oc -A 



2kl + k 2 y 



16 



t-2 
n 



(Al) 



This implies that the tree-level gap is actually filled. Does 
this invalidate the step of integrating out the field if as 
done in Sec. 1111 Al ? In order to answer that question one 
has to estimate the emerging interaction terms. The most 
relevant one of those is in the static limit given by 



SS< 



d d xdr 



with 



d d qdu! 



0,w n = 0) (A2) 



Xffff( fc = 0,w„ = 0) oc J / 27[ N d+1 Gn(q> u)G„(-q, -u>) 

(A3) 

in which G^(g, w) now is the full Green's function taking 
into account the self-energy in Eq. (|A1|) . The respective 
integral reads 



X^{k = 0,u) n = 0) oc 



d D qduj 



1 



(27T) 



(A4) 



with a parameterizing the strength of the sclf-cncrgy cor- 
rection. This expression reduces to case of gapped if for 
a = 0. It is obvious that a^O docs not produce singu- 
lar contributions. We have performed similar checks for 
other terms generated from integrating out if. 

We thus conclude that the presence of the cubic term 
(fT4|) . although rendering the field 7? gapless at critical- 
ity, does not invalidate the derivation of the effective co 4 
theory ([IS]). 



Appendix B: Derivation of the non-linear sigma 
model 



Here we sketch the derivation of a non-linear sigma 
model for the staggered-dimcr Hcisenbcrg antifcrromag- 
net, as usual performed in the scmiclassical limit start- 
ing from the magnetically ordered phase (see Chap. 13 
of Ref . 0) . A cubic term will appear as a result of this 
derivation, further discussed in Sec. MI Bl 

The Hamiltonian ([3]) can be written as 



H — ]P J,aS, ■ Sj+s, 
38 



(Bl) 



where 5 = x,y and 

J j5 = J [1 + (A/2)&, 4 (1 + (-1)'')] . (B2) 



Here, A = J'/ J — 1 measures the modulation in the 
couplings, and (— 1) J = ±1 for the two sublattices of the 
square lattice (solid and open circles in Fig. [I]}). 

To derive a field theory, we replace Sj —> SNj, where 
Nj is a three-component unit-length vector. Assuming 
proximity to a state with collinear Ncel order, Nj can be 
parameterized by 



Ni(r) = {-iyftj{r) 



1 - &%<r) 



1/2 



(B3) 

Here and Li are the (slowly varying) staggered and 
uniform components of the magnetization, respectively, 
obeying the constraints rif = 1 and Hi ■ Li = 0. We have 
restored the lattice constant a and assume that Lf <C 
S 2 a~\ 

Substituting Eq. (|B3|) into the Hamiltonian (|Blj) and 
expanding the square root results in the following Hamil- 
tonian piece Jiy for the magnetic couplings along the y 
axis: 



2« 4 f 2 



L j+y 



L, 



+ (-1) 



3+Vf 



L a - 



0(L 3 



a 

(B4) 



Upon taking the continuum limit, only the first two terms 
of Eq. (|B4[) are finite. The other two terms oscillate on 
the lattice scale and disappear in the continuum limit. 
The remaining Hamiltonian piece, H x , is similar to T-L y 
with y — > x, with the crucial difference that the oscillat- 
ing behavior of the couplings changes the prefactors of 
the third and fourth term to (— 1) J (— l) 3 = 1. Therefore, 
these terms - which will eventually lead to a cubic term 
analogous to Eq. - survive in the continuum limit 
for the staggered-dimer model (but not for the colum- 
nar dimer model, simply reflecting that a cubic term is 
forbidden by momentum conservation in the latter). 

The continuum version of the Hamiltonian (|B1[) then 
reads 



H 



d 2 rdr {S 2 [(d y n) 2 + (1 + A/2) {d x n) 2 } 

+ 4a 2 (2 + A/2)L 2 -2SaAL- (9 x n)| (B5) 

for the staggered-dimcr model, while the last term is ab- 
sent in the columnar dimer model. 

Passing to a coherent-state path-integral formulation 
and integrating out the L fields, the action for the unit- 
length n field assumes the form 



where 



S — Sb + S2 + S3, 



,/3 ,1 

Sb = iS / ](— 1) J / dr I dutij ■ (d u Hj x d T tij) (B7) 
^0 Jo 



(B6) 
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is the familiar Berry phase term [Eq. (13.52) of Ref. Q 
and 



S 2 

S, 
with 



1 

. JSAa 



9 



d 2 rdr [c 2 x (d x n) 2 + c 2 y (d y n) 2 
d 2 rdr n ■ (d x n x d T n), 



(d T n) 



(B8) 



JSay/8 + 6A, 



JSaVs + 2A, 
2a 2 J(4 + A). 



The velocities c x and c y agree with the spin-wave ve- 
locities calculated within spin-wave theory for the model 
(1B1[) . The prcfactor of 53 vanishes for A = 0, i.e., for the 
(unmodulated) square-lattice Heiscnbcrg model. 



Appendix C: RG analysis of the non-linear sigma 
model 



Despite the fact that the cubic term (|B8[) in the non- 
linear sigma model involves an artificial quantization as 
discussed in Sec. 1111 B[ one can attempt a RG analysis 
of this non-linear sigma model using an expansion in e = 
(d— 1). Here, we follow the calculation of Rcfs. l46ll47l and 
only display the relevant changes. 

We assume an isotropic velocity c and a momentum- 
space ultraviolet cutoff A. After re-scaling the coordi- 
nates according to xq = Act and x' = Ax, the action 
assumes the form 



2<?o 



dxo / d d x [(d^fi) 2 — 2hgon z 



2i5gori ■ (don x din)] 



(CI) 



where u = 0,1,..., d, u — Ac/3 is the re-scaled tem- 
perature, l/go = p5/(cA d ~ 1 ) measures the stiffness ps, 
h = H / (cA^ d+1 ^ encodes an applied staggered field H, 
and 8 = A/(cA d_1 ) represents the strength of the cubic 
term. The tree-level scaling dimensions of the coupling 
constants follow as [go] = 1—d, [h] = d+l, and [8] = d—1. 
In contrast to the -theory, here power-counting indi- 
cates that the cubic term is relevant for d > 1. 

The e expansion is generated as usual via the parame- 
terization ft = (tt x , n y , [1 — n 2 — it 2 } 1 / 2 ) and expansion of 
the action in 7? fields. Doing so, the lowest-order contri- 
bution of the original cubic term is found to be of order 
0(tt 4 ). Momentum-shell RG equations are obtained from 



integrating out modes with momenta e~ l < k < 1 and 
diagrammatically analyzing the perturbative corrections 
as in Ref. It turns out that the new coupling S does 
not modify the one-loop flow of g and h, as the possible 
contributions exactly cancel. The only correction to 8 
arises from a 5g diagram. 

As a result, the one- loop RG equations in the limit 
H ->• and T -> read: 



dl = {d 



1)6- 



d+l 



KcigS, 



(C2) 
(C3) 



where A7 1 = 2 d " V/ 2 r(d/2). Eq. JC2]) corresponds to 
the limit T ->• of the equation (3.1a) from Ref. \It\. 

The renormalization group flow is illustrated in Fig. [7] 
The only non-trivial fixed point is at (g, 8) = (g c , 0) with 
g c = 2{d-l)/K d . While this controls the QPT for 6 = 0, 
it is unstable w.r.t. finite 6. While it is possible that 
the inclusion of higher loop orders stabilizes a non-trivial 
fixed point at finite 6, the one-loop result in itself is puz- 
zling: Most disturbingly, the coupling 6 becomes more 
relevant with increasing d (already at tree level), in con- 
trast to conventional expectations. One might argue that 
in fact the combination (g6 ) , being marginal at tree level, 
measures the strength of the cubic term. However, the 
absence of a stable fixed point remains to be understood. 

Thus, the RG results for the non-linear sigma model 
(above) and for the ip 4 model (Sec. I IV A[) appear to mu- 
tually disagree regarding the role of the two-particle de- 
cay term. As discussed in Sec. IIII B[ we believe that the 
non-linear sigma model analysis is not trustworthy, (at 
least partially) related artifacts of the unit-length con- 
tinuum limit. It is worth mentioning that disagreement 
between the two field theories was already pointed out, 
for instance, in Refs. liH andl49l. To our knowledge, these 
issues are not completely settled^ 




FIG. 7: Schematic renormalization group flow for the nonlin- 
ear a model (|C1[) for the case d > 1, h = 0, and T = 0. 
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